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Abstract 

We derive a set of equations expressing the parameters of the magnetic inter¬ 
actions characterizing a strongly correlated electronic system in terms of single¬ 
electron Green’s functions and self-energies. This allows to establish a mapping 
between the initial electronic system and a spin model including up to quadratic 
interactions between the effective spins, with a general interaction (exchange) 
tensor that accounts for anisotropic exchange, Dzyaloshinskii-Moriya interac¬ 
tion and other symmetric terms such as dipole-dipole interaction. We present 
the formulas in a format that can be used for computations via Dynamical Mean 
Field Theory algorithms. 
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1. Introduction 


Describing a solid in terms of its magnetic properties requires the knowledge 
of an effective spin model which displays the same interesting physical properties 
as the many-electron Hamiltonian whose exact solution would give the complete 
description of the system. The determination of the form of the effective spin 
model and of the strength of the interactions between the constituent spins 
starting from the initial electronic model is, in general, a complicated many- 

body problem [laiiiiSiiEEllIl. 

We have recently derived expressions for the parameters of the magnetic in¬ 
teractions within an extended (multi-orbital) Hubbard model , in the pres¬ 
ence of arbitrary relativistic couplings affecting the electronic degrees of freedom 


* Corresponding author 

Email addresses: a.secchi@science.ru.nl (A. Secchi ), andrea.secchi@gmail.coin (A. 
Secchi ), alichten@physnet.uni-hamburg.de (A. 1. Lichtenstein), 
m.katsnelson@science.ru.nl (M. I. Katsnelson) 


Preprint submitted to Elsevier 


June 23, 2015 








(such as spin-orbit, magnetic anisotropy, Zeeman coupling with an external mag¬ 
netic field). The formulas presented in Ref.[I^, after neglecting the vertices of 
two-electron Green’s functions, are expressed in terms of single-electron (but 
fully interacting) Green’s functions G and the sin gle- electron (hopping) Hamil¬ 
tonian T. The use of the representation via T [9|, [l^ for computations related 
to real materials requires the additional step of a tight-binding parametrization, 
which is implemented only in some methods of electronic structure calculations. 
On the other hand, a presentation of the formulas in terms of Green’s functions 
G and self-energies S would make them more suitable for implementation via 
Dynamical Mean Field Theory (DMFT) [l^, 14, 15|, since any DMFT calcu¬ 
lation deals with G and E. Writing the parameters in a ww that explicitly 
exhibits self-energies, analogous to what was done in Refs.|4 Isl. Il0l|. also allows 
to explicitly include the approximation of local self-energy, which is the key 
assumption of DMFT. We here present the adaptation of the formulas for the 
exchange tensor to this scheme. 

2. Method and discussion 


We consider the extended multi-orbital Hubbard Hamiltonian [16|, [17 

ilHQ. 
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where the field operator ^ m creates an electron with quantum numbers 
{o, cr, m}: o refers to a set of the orbital indices (for a basis of localized Wan- 
nier wave functions, these are the atom index a, the principal atomic quantum 
number n and the angular momentum quantum number 1: o = {a,n,l}), while 
cr € {t, i} and m € {—I, + are the third components of the intrinsic- 

spin and orbital angular momenta, respectively. Local angular momenta are 
measured with respect to local reference frames, which depend on o and might 
not be collinear [l2|. The single-particle Hamiltonian matrix is com¬ 

pletely arbitrary, so it can include any relativistic single-electron terms (Zeeman 
coupling, spin-orbit, magnetic anisotro pies ). The interaction Hamiltonian Hy 


is assumed to be rotationally invariant 12|. 


The goal in Ref. 12| was to map the model given by Eq. ([1]) onto an effective 
model of classical spins Bq including up to (arbitrary) quadratic interactions, 
with Hamiltonian 
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determined by the exchange tensor 77“^ = 77“ “ (here and in the following 
a and a' are used to denote the space coordinates, e.g. x, ?/, z) and the ef¬ 
fective magnetic field l&o- It is convenient to decompose the exchange ten¬ 
sor into the three vectors 77oo' = o'o (anisotropic exchange), 'Doo' = —T^o'o 
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(Dzyaloshinskii-Moriya interaction), and Coo' = Co'o (symmetric non-diagonal 
exchange), defined as 
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where e““ “ is the completely anti-symmetric tensor of rank 3. The Heisenberg 
model is obtained as the particular case in which Joo' ■ 

To perform the mapping, in Ref. 1^ we have derived the response of the ther¬ 
modynamic potential of the electronic system under small spatially-dependent 
rotations of the spin quantization axes associated with each orbital spinor de¬ 
noted by o, up to second order in the rotation angles. The derivation of such 
response involves path integration over the fermionic fields after the introduc¬ 
tion of auxiliary bosonic degrees of freedom which express the amplitudes of 
rotations from an initial spin configuration; the coefficients of the interactions 
between the remaining bosons are put in correspondence with the parameters 
of the spin model Q by imposing that the thermodynamic potential of the spin 
system after the spin rotations is equal to that of the electrons. Excluding the 
vertex contributions, the parameters of the spin model are expressed in terms 
of single-electron Green’s functions (which of course include interaction effects) 
and the single-particle part of the electronic Hamiltonian, T. 

This procedure is similar to the one previously adopted in Refs. ii for the 
case of quenched orbital moments, but in Ref .[13 we have considered rotations 
of the local total spins So = lo + So, where lo and So are, respectively, the orbital 
and intrinsic angular momenta associated with the states o. More precisely, we 
have considered rotations in the space of the single-particle eigenfunctions of 
So and So, analogously to Ref.Q, while in Refs. [1,1^ the rotations affected 
the space of eigenfunctions of and s^. This allowed us to obtain formulas 
for the exchange tensor that can be separated into contributions coming from 
the interactions between spin-spin, orbital-orbital, or spin-orbital degrees of 
freedom of the electrons. It should be noted that this possibility is not applicable 
within Density Functional Theory (DFT) formulations, where observables are 
expressed in terms of the charge density and the intrinsic-spin density. The 
possibility of rotating local total spins is related to the representation of the 
electronic Hamiltonian in terms of localized wave functions, which implies a 
higher number of degrees of freedom with respect to DFT (related to the fact 
that the set of localized states would be over-complete in theory, or not even 
complete in practice due to truncation). 

The computation of the magnetic parameters via DMFT is greatly simplified 
if they are formulated in terms of single-particle Green’s functions and self¬ 
energies F in magnetically ordered states, since this avoids the initial step of a 
tight-binding parameterization of the single-electron Hamiltonian T. To remove 
T and introduce E, we use the equations of motion for Matsubara Green’s 
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functions (Dyson equations), which we write in general matrix notation as 

(w — i/r) G(iw) + iT ■ G(iw) = 1 — S(iw) • G(ia;), 

{lo — i/r) G(ia;) + iG(iw) • T = 1 — G(ia;) • E(ia;). (4) 


These equations hold for the Matsubara Green’s functions defined according to 
the following convention: 

Glir) ^ -i(r^i(r)^|) ^ i^Gl(ia;)e--^ ( 5 ) 

^ (jJ 

where w = (2n + l)7r//3 is a fermionic Matsubara frequency. As a particular 
case, the single-electron density matrix is given by 

P =-iG(r = 0 ”) =-i^^e“^°'*'G(ia;). ( 6 ) 


We now have to distinguish between the magnetic parameters that can be 
computed from the second-order response in the rotation ang les and those which 
are computed from the first-order response. From Ref.[l2l|. we note that the 
former terms can all be written in terms of the following quantity: 
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where X, Y G {spin, orb} refer to either spin- or orbital- related terms, that is. 
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where Soa is an intrinsic spin matrix (uoa is a Pauli matrix), while loa is 
an orbital angular momentum matrix. In Eq.([7]) we have used the notation 
{A; B} = A ■ B + B ■ A to denote the anti-commutator of the matrices A and 
B\ in the following we will also make use oi\A\B] = A - B — B ■ A to denote the 
commutator. 

From the Dyson equations ([4]), we have (the frequency arguments of Green’s 
functions G and self-energies E are implicit): 


T • G = -i [1 - (w - i/r) G - E • G], G ■ T = -i [1 - (w - i/r) G - G ■ E], 

T • G • T = -iT -h E - E • G • E -h (w - i/r) (1 - E • G - G • E) - (w - i^)^ G, 
[T;p]=Tr^[E;G], (9) 
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where we have introduced the notation 

1 = Tr^(j .g). 

OJ 

Applying Eqs.® to Eq.(I7]), we obtain 
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W • {E;G}:) 
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- Si, • G°, • sj,,,, • [E. G • qf} 
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We then consider the magnetic parameters determined from the first-order 
response. From Eqs.(68) of Ref. 12|, we see that these are m 

Cqj,,, and Cl,. The first-order response term (in the RHS of Eqs.(68) of Ref. [3) 
can be written as 


= iTr™,, {Si • [p;r]:;) = iTr^.<,Tr^ {Sl • [G-E]°l 
= iTr^,.Tr^ ^ • (g°, ■ E°' - E°, ■ G°')' . 

o' 


( 12 ) 


By separating local and non-local terms, as well as taking into account the 
symmetries of the latter, it is then possible to identify the remaining parameters 
of the spin model. It should be noted that the parameters obtained with this 
procedure are not equivalent to those expressed in terms of the single-electron 
Hamiltonian T in Refs. 0 and Q. These are different definitions, which respect 
the defining equations (68) of Ref.[l3, but are more directly applicable for a 
DMFT implementation. 

In the next section we list the resulting formulas for the magnetic parameters. 


3. Results 

3.1. Dzyaloshinskii-Moriya interaction 

The Dzyaloshinskii-Moriya parameters, given by the vector 'Dqo' > are written 
as 


{Vir'^ = ^Tr^,^Tr^ • (g°, • E°' - E°, • G°') 


(V^ — -Tr Tr 
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= ^Tr„..Tr^ [soy ' (g°, • S°' - E°, • G°') 

- So,, • (g°' • S°, - E°' . G°,) ■ , 

= ^Tr^.,.Tr^ [loy ' (g°, • E°' - S°, • G°') 

-lo-y(Gi -K'-K -Gl)]. (14) 

^jyz jspin-spin 

2 \orb—orb 
oo') 

fj^z jspin-orb 

3.2. Symmetric out-of-diagonal interactions 

The symmetric out-of-diagonal interaction parameters, given by the vector 
Coo' with o ^ o', are written as 

{Cfo'r^^ = ^Tr™,.Tr^ [so. • (g°, • E°' - E°, • G°') 

+ So,, • (g°' • E°, - S°' • G°,) ], 

Ko'T'^ = ^Tr^.^Tro, [lo. ■ (g°, • E°' - E°, • G°') 

+ /o'.- (g^E°,-E^G°,)], (16) 
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(0)^"“ = -^Tr™..Tr<o [so, • (g°, • E°' - E°, ■ G°') 

+ So,,-(g^E°,-E°'.G°,)], 

(C')”" = -^Tr„.,.Tr^ ^oy ■ (g°, • E°' - E°, • G°') 

+ Zo',- (g^E°,-E°'.G°,)], (17) 
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3.3. Local out-of-diagonal anisotropy 

The local out-of-diagonal anisotropy is given by the vector Coo- We have 


\spin —spin _ TTispin, spin 
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(^.^)spi„ ^ iTr^Tr^_^ 1 ]°]) - , 

{CfoT'^ = i Tr<,Tr„,. • [G°; E°]) - , (20) 


(^,^)spi„ ^ _iTr^Tr„_^ . [gO; E°]) _ , 

(a)”" = -iTr^Tr^,. {Ly • [G°; S°]) - {B^T^ , (21) 

where the components of the effective magnetic field are 

" = 5i/2A^bBo [Tr, (soTr„p°) • <], 

= gifiBBo [Tr„ (bTr,p°) • <], (22) 

with gi /2 and gi being the intrinsic-spin and orbital (/-factors, respectively, and 
Bo the value of the external magnetic field acting at the position of the orbitals 

o. 


3 . 4 . Anisotropic exchange interactions 

The anisotropic exchange parameters are given by the vector 4J 00 ' with 
o o'. They are written as 
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and the terms related to are given by the averages of the respective terms 

contributing to and since /2. 

It is instructive to consider the particular case where the self-energy is local 
not only in (atom) position space, but also diagonal in the principal quan¬ 
tum number (n) and orbital angular momentum (1) indices (we recall that 
o = {a, n, 1}), so that 


S°, « Sf,^„ 
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we obtain 



—spin ^ 

i -Tr 

cnTr,.^ 

,a [^oy 5 ^o] 

•G°, 





■orb ^ 

-Tr<^ 

Tlm,CT 

' ^ [^oy: ^o] 

/^O 

^ o' 

\lo'y: Sq/] 

g:}, 


{JSo^ 

—orb ^ 

_T~'t T~'t 

-Licj 

y ['^oy; ^o] 

•G°, 

[^o^y 5 

•G°' 




+ [So 

'y^ ^o' 

’]-G^'-[loy 

;So]' 



(26) 


-spin 

^ -Tr, 

Tr 

uj -Lim, 

,<7 [^OX -i ^o] 


[^o'a: 5 ^o' 




•orb ^ 

-Tr<^ 

Tr 

-Li m,o' 

1 \}‘Ox] ^o] * 

G^,- 

^o'] * 

G°'}, 



-orb ^ 

_T't' T~'t 

-Liu; -Lim,{ 

7 ^ \_^OX t ^o] 

•G°, 

U'o'x ; 

•G°' 




+ [So 

'x 5 ^o' 

’] • • [^Otc 

;So] 

■ g:,}. 


(27) 

ighly symmetric 

case, 

also assuming collinear magnetic states and that 


Green’s functions and self-energies are diagonal in spin space, one obtains the 
isotropic exchange parameter 
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as 


Tr^Tr™ (sf • • Sf, • + Sf • G°f^ • Sf, • G^^) , (28) 


where Sf = (l]| — E);) /2, which is consistent with Eq.(109) of Ref. 12| and 
with the previous literature SS0- 


3.5. Local exchange interaetions (diagonal anisotropy) 

The local exchange interaction parameters, or diagonal components of the 
local magnetic anisotropy, are given by Soo- Within our rotational procedure, 
it is possible to determine only two of the three parameters as a function of the 
third one. Putting (a, a) = {x,y) or {y,x), we obtain 

/ /T-a \spin—spin j^spin, spin i /'ioj2\spin . f \ ' / zzx . zj-y \spin—spin 

oo ^oo) ^ooL^oa. ' V^o/ ' 2 / x 00' ^00't / t 

o' 

t q-OL _ qz \ orb—orb _ porb,orb 1 ('R^'vorb . ^ \ '' / qx 1 qV \orb—orb 

00 ^00) ^oci^oci ' \^o) 2 / ^ 00' ' ^00') "> 
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/ q^ _ qz \ spin—orb _ ^spin, orb , ^orb, spin , ^ \ / qx i qV \spin—orb 

^00) '^oct,oct ^ooL^ooL ' / j oo' ' ^oo') ’ 

o' ^o 

(29) 

where the non-local exchange terms are given in Section [3.41 and the magnetic 
field components are defined in Eos. ([2^ . 
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4. Conclusions 


We have provided the formulas for the general exchange tensor expressing 
the quadratic magnetic interactions in strongly correlated systems, in a version 
that can be implemented via DMFT. The formulas allow to compute the effects 
due to intrinsic-spin and orbital degrees of freedom of the electrons on equal 
footing (the orbital magnetic moments are not quenched), and to distinguish 
between spin, orbital and spin-orbital interactions that contribute to the ex¬ 
change tensor. The obtained formulas represent the extension to the relativistic 
case and the generalization to unquenched orbital mametic moments of the well- 
known formulas for spin-only exchange interactions [^, which are recovered as 
a particular case. We remark that effects due to the non-locality of self-energies 
in position space are included in our theory both as presented in Ref.[l^ and 
as presented here; although they cannot be computed within DMFT, a possible 
approach to include them is via the Dual-Fermion scheme 


211 . 
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